Dissipative dynamics of spins in quantum dots 
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We present a theory for the dissipation of electronic spins trapped in quantum dots due to their 
coupling to the host lattice acoustic phonon modes. Based on the theory of dissipative two level 
systems for the spin dynamics, we derive a relation between the spin dissipative bath, the electron 
confinement, and the electron-phonon interaction. We find that there is an energy scale, typically 
smaller than the electronic lateral confinement energy, which sets the boundary between different 
dissipative regimes . 



Manipulating quantum states of matter to achieve new 
ways of processing information has been a subject of in- 
tense research in the last decade 0. There is now a 
wide variety of proposals to take advantage of the inter- 
nal transformations of quantum systems to perform the 
so-called quantum information processing (QIP) Q. In 
particular, condensed matter systems such as SQUIDS 
and semiconductor quantum dots are seen as possible 
candidates for future implementations of QIP devices. 
In these systems, it is possible to find or create well de- 
fined two-level quantum states that potentially fulfill the 
minima criteria proposed by Di Vicenzo to be quantum 
bits (qubits) 0. However, apart from the technical dif- 
ficulties to implement acceptable qubits, nature imposes 
another limitation which is ubiquitous. This comes from 
the impossibility of perfectly isolating a complex quan- 
tum system from its environment that results in the loss 
of quantum coherence 0, @|. 

In this paper we investigate the decoherence of elec- 
tronic spins in quantum dots due to their coupling to 
the acoustic phonon modes of the host lattice. In few- 
electron quantum dots the electronic spin is a good can- 
didate for a qubit and the confinement of the electronic 
wave function plays a major role in isolating the spins 
from most energy relaxation channels. In fact, using 
a second order perturbation approach, Khaetskii and 
Nazarov 0| have evaluated the spin-fiip transition rate 
between Zeeman sublevels due to its coupling to acous- 
tic phonons. They showed that, since angular momen- 
tum conservation requires the spin-fiip process to involve 
a virtual transition between excited orbital states, the 
spin-fiip rates are suppressed as ujq'^, where ujo is the lat- 
eral confinement energy. According to their result, even 
for Wo ~ 1 meV , the spin-fiip rates can be of the order 
of milliseconds. 

Quantum dots can be fabricated by confining a two 
dimensional electron gas of a semiconductor heterostruc- 
ture into a region of the order of the Fermi wavelength. 
In the "lateral quantum dots" (LQD) this confinement is 
done laterally, by means of surface electrostatic metallic 
gates. In "vertical quantum dots" (VQD) the electron 



gas is confined vertically by etching techniques that cre- 
ate a circular pillar heterostructure. The typical lateral 
confinement length achieved with these fabrication tech- 
niques are of the order of hundreds of nanometers, which 
leads to wqw 1 meV for LQD Q and ~ 3 — 5 meV for 
the VQD Much higher confinement energies can be 
achieved in self-assembled quantum dots (SAQD) which 
consist of ensembles of dislocation free semiconductor 
nano-crystals embedded into a semiconducting matrix of 
different band gap For InAs:GaAs SAQD ^ the 
electronic wave function extent is of the order of 50 A, 
resulting in loq ~ 50 meV. For all the above systems 
single electron charging is easily observable due to their 
large Coulomb blockade gap, which can be of the order 
of 2 meV for LQD and VQD and as large as 20 meV for 
SAQD. 

Although no experiment has so far been able to deter- 
mine directly the decoherence rate, of spins in quan- 
tum dots, there are several attempts to indirectly deter- 
mine lower bounds to the spin relaxation rate Pail- 
lard and co-workers have studied the time-resolved 
photoluminescence of InAsiGaAs SAQD and observed 
that, within the time scale of an exciton lifetime, the 
carrier spins are totally frozen. Fujisawa and co-workers 
[1| have demonstrated that orbital sublevel transitions of 
VQD involving spin-fiips have relaxation times that are 
4 to 5 orders of magnitude longer than those which do 
not involve spin-fiips. This gives an idea of the degree of 
isolation of the spins in the VQD. More recently, Hanson 
and co-workers [3| used short voltage pulse sequences to 
measure the relaxation time of spins in LQD and con- 
cluded that the lower bound of Ti is 50^s for a magnetic 
field of 7.5 T. 

In quantum dots, the most important channels of dis- 
sipation for the spin are indirect and via the spin-orbit 
interaction. The orbital dissipative dynamics on its turn 
is dominated by the electron-phonon interactions. It is 
thus desirable to have a theoretical description which can 
connect the spin and orbital dissipation channels. A sim- 
ilar problem of indirect dissipation appears on a com- 
pletely different context of electron tunneling between 
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the atoms of a diatomic molecule embedded in a viscous 
environment Here we analyze the spin dissipation 

in a quantum dot from the same perspective. We de- 
rive an effective "bath" spectral density seen by the spins 
that results from the spin-orbit coupling and the orbital 
damping. 

Following previous theoretical works we treat the elec- 
trons in the effective mass approximation and consider 
the confining potential for the envelope wave function as 
paraboHc. This was shown to be a good approxi- 
mation for the LQD and VQD with low electronic fill- 
ing. Also, for SAQD this phenomenological model de- 
scribes the orbital electronic density of states probed by 
magneto-capacitance measurements with fine accuracy 
up to the 3th excited level of the dot We further 

assume that the harmonic frequency in the direction per- 
pendicular to the quantum dot plane, ll!±, is much higher 
than the lateral harmonic frequency ujq. Therefore the 
relevant low energy orbital degrees of freedom are in the 
X — y plane whereas the orbital dynamics in the z direc- 
tion is practically frozen. This is a reasonable assumption 
even for typical InAsiGaAs SAQD since ^ « 8 - 10 
For the spin degrees of freedom, besides the Zee- 
man term arising from an external magnetic field in the 
z direction, we include the Dresselhaus spin-orbit interac- 
tion projected on the x, y plane 0| which is responsible 
for the coupling between the spins and the dissipative 
phonon bath. Thus, apart from a zero point energy on 
the z direction, the spin-orbit Hamiltonian will be 
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where A = giJ-sBz (we use h = 1), (3 = ^rj^l) = 
"fc'm*uj_L (with 7c being the Kane parameter m* 
is the electron effective mass and g its gyromagnetic fac- 
tor. The operators a^^y) are the usual ladder operator 
for the X (y) direction. Notice that although we set an 
external field for the spins, we are neglecting any dia- 
magnetic contribution for the orbital degrees of freedom 
like in the Fock-Darwin description. This simplification 
allows us to separate the degrees of freedom in the x 
and y directions and can be well justified for ^ 
where ujc = eB/ {m*c). One should also notice that, as 
opposed to the notation on reference 0, A here plays 
the role of a "tunnehng" field rather than the "bias" field, 
even though it refers to the z direction. This happens 
because dissipation occurs only in the x — y plane of the 
dot. 

The electron-phonon coupling in this restricted sub- 
space can be written as 
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where C^ x is the electron-phonon coupling for phonons 
with polarization A and frequency > , and r = {x, y, z) 
is the electron position operator. Here we consider only 
the piezoelectric and deformation potential interactions 
with acoustic phonon modes in zinc-blende structures 
It can be shown that within the linear response ap- 
proximation for the phonon system the electron-phonon 
Hamiltonian (EqlJl is mapped into the hath of oscillators 
model with the spectral function given by 



Js (w) 



where s = 3 for the piezoelectric interaction, with dimen- 
sionless coupling 5^ = ^^"^I'l"-" ( ^ _|_ J^. ) a^d s = 5 



for the deformation potential, with ^5 



where 



cj£) is the Debye frequency, vi and Vt are the longitudi- 
nal and transverse sound velocities respectively, p is the 
material density, (em) 14 is the electromechanical tensor 
for zinc-blende structures Ii.tI. a nd Oct is the deforma- 
tion potential in the F point is the Heaviside step 
function. 

Now, since the spin degree of freedom is coupled to the 
orbital motion of the electron, we can adopt the prescrip- 
tion of reference to extract the spectral function of 
the effective heat bath to which the spin is now coupled. 
As opposed to the case of reference where both the 
"spin" and the phonon degrees of freedom are coupled to 
the position of the orbit, here we have the phonon coupled 
to the position (Eq[2l and the spin to the momentum of 
the electron (EqQ . This leads to a significant change on 
the "effective bath" spectral function, namely 



s+2 



' {cod - to) , (3) 



where Z {^) ^ (^) - (^) (l + 5, 0, (-i^) ) , 

and (j)s (x) = \V .t"/ (y^ - yx'^) dy = 
-^{B{x,s,0) + {-iy B{~x,s,0)), with B be- 
ing the generalized incomplete beta function. For x <C 1 
we can approximate 0s (x) ~ ^ (7^ + 

For most semiconductors used in the fabrication of 
quantum dots the dimensionless constants 63 and 65 are 
of the order of 10^ and 10^ respectively. For instance, 
using the bulk physical paramenters we obtain, in 
GaAs S3 = 355 and S5 = 1.95 x 10*^ and in InAs 
S3 = 149 and S5 = 5.03 x 10'^. Furthermore, the typi- 
cal frequency for the spin dynamics is much smaller than 
cj_D « 30 — 50mey, which suggests an asymptotic analy- 
sis of Jeff (uj). In order to do that we should first notice 
that Jeff (w), Eq. |S1 is peaked at uj = fls, where fls is 
defined as the solution of Z (fls) = 0. For ujq < lod this 
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solution is given by: 
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In the weak coupling limit, i.e., Ss <C 1, ils ~ 
(^1 — ^(^L2) ) 3,nd the resonance at fig corresponds to 
the harmonic frequency of the orbital dynamics, shifted 
by the electron-phonon coupling. However, for (5s ^ 1, 
the peak originally centered at the lateral confinement 

frequency, ujq, is drastically shifted to Qs =^^01 1 ^'^^'^ 
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In the low frequency range, defined by w <C 
— < — ^ , we obtain 



and 
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Note that in this limit the effective spectral function felt 
by the spins is always super-ohmic, with a power s + 2. 
On the other hand, in the high frequency limit, Vis ^ 
<C the spectral function can be approximated by 
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In this limit the spectral function has a power s — 2 
that can even be sub-ohmic. The behavior of Jeff ((^) 
is sketched in figure ^ for the piezoelectric interaction 
(. = 3). 

All this complex structure of the effective spectral func- 
tion will provide us with a new physics for the system. 
Within the adiabatic renormalization scheme, the fast 
modes of the bath, i.e. those with frequency much higher 
than A, can quickly adjust to the spin fiip motion and 
are successively integrated out by a Born-Oppenheimer 
approximation 

At long times the spin fiips coherently with a renormal- 
ized Zeeman frequency A^, dressed by a Franck-Condon 
(FC) factor, which represents the polarization cloud of 
the high frequency phonon-orbit complex, given by : 
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where F {z) = ^ duj^^^^j^ coth (^) is the FC factor 
with p being an unimportant dimensionless number 
much larger than 1. In the low frequency limit z ^ 
fls the FC factor is dominated by the region u fts 
where '^^ff,j'^^ can be well approximated by a Lorentzian 
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and thus 
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. Since A^ < 



A, the solution for A ^ fig involves integration over the 
entire spectrum of Jeff {(jj), where most of the weight is 
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Figure 1: Full lines, exact Je// {uj) (s = 3, ^f- = 0.5 and 
^3 — 100), dashed lines, asymptotic behavior according to 
Eq. IHand Eq. H (inset). 



concentrated around the peak centered at fis, and thus 
the lower limit of integration can be extended to zero, 
yielding 
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For temperatures T ^ we then 



have Ar ~ A exp | — ^ | . From this point onwards we re- 
strict our discussion to the piezoelectric electron-phonon 
coupling which leads to a larger bath spectral density at 
low frequencies. In the low frequency limit, defined by 
A ^ r^a , the dynamics is dominated by a super-ohmic 
relaxation with a power s + 2 (see Eq. 31- This allows 
the spin to present coherent damped oscillations with a 



decoherence rate 
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Except for the renormalized Zeeman frequency A^, 
this is essentially the spin-fiip rate from Khaetskii and 
Nazarov [6] which, for GaAs quantum dots with m*/3^ ^ 
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ImeV and A = 0.025 meV, 



gives decoherence times of the order of milliseconds . One 
should note however that, since fl^ ^ O (0.1 wo), for Zee- 
man splittings in the range A ^ 0.1 — 1 meV the pertur- 
bative result is only applicable for dots with lateral con- 
finement energies uq ^ 1 meV. This is usually the case 
then for SAQD, in which one can have cuq ~ 50meV. 
For the LQD and VQD wq ~ O (1 meV) and thus A can 
go beyond this perturbative limit. For A ~ f^a, the re- 
laxation is dominated by the resonance of the spin with 
the orbit-plus-phonons composite at frequency fls. This 
resonance has a linewidth of £3 ~ (^) W' which qual- 
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Figure 2: Density map of Jeff / (rn* (3^^ from piezoelectric 
interaction (s = 3) for GaAs (left) and InAs (right) QD. 

itatively determines the decoherence rate of the spin in 
this intermediate range of frequencies. 

In the higher frequency range fig ^ A ^ cj_d the inte- 
gral on equation © is dominated by the high frequency 
part of Jeff (uj). For s > 3 this integral is algebraic and 
only leads to a small correction. However, if s < 3 the 
behavior is qualitatively different and dominated by the 
lower limit of the integration which has an "infrared di- 
vergence". For the piezoelectric interaction (s ~ 3), it 
behaves at low temperatures T <C Od as 

where K = "^^ ^ ^g^^^ <C 1 . In this "high frequency 
region", fi^ ^ A ^ wd, that can be achieved for wo ^ 
ImeV ("large" quantum dots), the decoherence will be 
dominated by an ohmic-like dissipation given by Eq. 
which yields a decoherence rate 

T2 26s \i^D ) 

Figure [2 shows the behavior of Jg/ / for GaAs and 
InAs quantum dots. The darkest region to the left 
(right) of the brighter stripe (corresponding to the peak) 
is the super-ohmic (ohmic) dissipation. A perpendicular 
magnetic field will increase the lateral confinement ac- 
cording to the Fock-Darwin energy ojfd = + I'^c- 
This is the case for the experiment on reference p] where 
the magnetic fields are of the order of 10 T and thus the 
cyclotron frequency ujc contribution can be higher than 
LUo- The density map of Figure El is then a guide for one 
to reach lowest decoherence rates given the lateral con- 
finement energy of the dot and the Zeeman splitting A. 

Thus, we conclude that the piezoelectric electron- 
phonon coupling only leads to large decoherence times, 
as predicted by the second-order perturbation theory, if 
the — ^ — i=. Otherwise the decoherence rate has a 
completely different behavior leading to much higher de- 



coherence rates. This suggests that as far as the decoher- 
ence due to acoustic phonons is concerned, SAQD better 
decouple the electronic spin degree of freedom from the 
environment. Nevertheless one should note that there are 
other decoherence mechanisms due to phonons. In ref- 
erence other electron phonon mechanisms were con- 
sidered extending the analysis of Khaetskii and Nazarov 
0. It has also been argued 113,11^01 that in SAQD, 
a resonance between the lateral confinement frequency 
and the longitudinal optical (LO) phonon mode can lead 
to a more efficient channel of orbital energy dissipation. 
In this scenario the orbital relaxation rate should be de- 
termined by the LO relaxation rate Fop (in bulk GaAs 
T~p « 7ps). Since the approach developed here con- 
nects the spin dissipation to the orbital dissipation and 
thus must also apply to these other orbital dissipation 
mechanisms if we replace J {uj) by the appropriate bath 
spectral function. 

This research was partially supported by Hewlett- 
Packard Brazil. The authors also acknowledge par- 
tial support from Conselho Nacional de Desenvolvimento 
Cientifico e Tecnologico (CNPq) and Fundagao de Am- 
paro a Pesquisa no Estado de Sao Paulo (FAPESP). 



[1] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information, (Cambridge University Press, 
New York, 2000) 

[2] D. Loss and D. P. Di Vicenzo, Phys. Rev. A 57, 120 
(1998) 

[3] D. P. Di Vicenzo, Science 270, 255 (1995) 

[4] A.O. Caldeira and A. J. Leggett, Annals of Physics 149, 

374 (1983) and Phys. Rev. A 31, 1059 (1985) 
[5] A. J. Leggett et al, Rev. Mod. Phys. 59 (1987) 
[6] A. V. Khaetskii and Y. V. Nazarov, Phys. Rev. B 61, 

12639 (2000), Phys. Rev. B 64, 125316 (2001) 
[7] R. Hanson et al, Phys. Rev. Lett. 91, 196802 (2003) 
[8] T. Fujisawa et al, Nature (London) 419, 278 (2002). 
[9] D. Bimberg, M. Grundmann, N. N. Ledentsov, Quantum 

Dot Heterostructures (.John Wiley & Sons, New York, 

1999) 

[10] M. Fricke et al Europhys. Lett, 36 (3), 197 (1996) 
[11] M. Paillard et al, Phys. Rev. Lett. 86, 1634 (2001). 
[12] A. Garg et al, J. Chem. Phys. 83, 4491 (1985) 
[13] P. Matagne et al Phys. Rev. B. 65, 085325 (2002) 
[14] E.L. Ivchenko and G. Pikus, Superlattices and Other Mi- 

crostructures, (Springer, Berlin, 1995). 
[15] P. Y. Yu and M. Cardona, Fundamentals of Semiconduc- 
tors, (Springer- Verlag, New York, 2001) 
[16] L. M. Woods et al, Phys. Rev. B 66, 161318 (2002) 
[17] X-Q. Li and Y. Arakawa, Phys. Rev. B 57, 12285 (1998), 
X-Q. Li, H. Nakayama, and Y. Arakawa, Phys. Rev. B 
59, 5069 (1999) 
[18] O. Verzelen et al, Phys. Rev. B 62, R4809 (2000) 
[19] S. Sauvage et al, Phys. Rev. Lett. 88, 177402 (2002) 



